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Abstract. We prove limit theorems for small-scale pair dispersion in synthetic velocity fields with 
power-law spatial spectra and wave-number dependent correlation times. These limit theorems 
are related to a family of generalized Richardson's laws with a limiting case corresponding to 
Richardson's t 3 and 4/3-laws. We also characterize a regime of positive dissipation of passive 
scalars. 

1. Introduction 

The celebrated Richardson's t 3 -law [SH] states that a pair of particles located at (x^ (t),x^ (t)) 6 
M. 2d being transported in the incompressible turbulence satisfies 

(1) E|x (1) (t) - x (0) (f)| 2 « C R et 3 for l x < - x (0) (t)| < £ 

where e is the energy dissipation rate, Cr is the Richardson constant and (.q and l\ are respectively 
the integral and viscous scales. Here and below E stands for the expectations w.r.t. the ensemble 
of the velocity fields. This law has been confirmed experimentally ( |22|. |39| . |31| ) and numerically 
( |43| . |18j . 0]). A stronger statement is that the relative diffusivity of the tracer particles 
is proportional to the 4/3 power of their momentary separation, and this is called Richardson's 
4/3-law (jSO]- See also (32], 0, EQ, [23 )■ This paper presents several small-scale limit theorems 
(Theorem m2 and 0) related to the Richardson's laws for a family of colored-noise-in-time velocity 
fields that have Kolmogorov-type spatial spectra and wave-number dependent correlation times. 
The other aspect of the scaling limit concerns the dissipation of the scalar field in the limit of 
vanishing molecular diffusion (Corollary 1 and 2). 

The nature of time correlation in fully developed turbulences in the inertial range is not entirely 
clear (see [3U] and the references therein). But it seems reasonable to assume that, to the leading 
order, the temporal correlation structure of the Eulerian velocity field u(t, x) is determined by the 
energy- containing velocity components above the integral scale, consistent with Taylor's hypothesis 
commonly used in the fluid flow measurements in the presence of a mean flow or the random 
sweeping hypothesis in the absence of a mean flow (see 00] and jSE])- In both cases the temporal 
correlation function on the small scales is anisotropic and depends on external forcing. The more 
robust features of small scale turbulence can be revealed by considering the relative velocity field 
U(t,x) = u(t,x + x(°'(t)) — u(t,x( ' (t)), with respect to a reference fluid particle x^°'(t), which 
tends to preserve invariance properties of the fluid equations. The velocity field u(t, x + x^°\t)) as 
viewed from a fluid particle, which is a useful tool for turbulence modeling [2], [212]; is called the 
quasi-Lagrangian velocity field in the physics literature and is an example of the general notion of 
the Lagrangian environment process [HE]) [HI], 



Department of Mathematics, University of California at Davis, Davis, CA 95616 Internet: fan- 
njian@math.ucdavis.edu . 



We assume ^3], ^H] that the two-time structure function of U(t,x) has the power-law form 
mU(t,x) - U(t,y)} ® [U(s,x) - U{s,y)\ 

2[l-cos(fc- {x -y))}e*v{-a\k\ w \t- s\)£ {hA)) {a,k)\k\ l ~ d dk, a 6(1,2), f3 > 0, a > 
with the energy spectrum 

(3) %, 4) (a,fc) - | Q) for |Jfc| ^ *o < oo, ^ > 0, E > 

where £\,£o are respectively the viscous and integral scales. The assumed temporally stationary 
vector field U(t, x) has homogeneous spatial increments and its expectation M s [U(t, x)}, conditioning 
on the events up to time s < t, is assumed to admit the spectral representation 

E s [U(t,x)-U(t,y)} 

(4) = / [l-exp(ifc- (x -y))]exp(-a\k\ 2 P\t - s\)U{s,dk), s<t 



where U (t, k) is a time-stationary process with uncorrelated increments over k such that 

(5) E[U(t,dk)U*(t,dk')] = £ {hM (a,k)6(k - k') dkdk', Vt,k,k'. 

The exponential form of the temporal correlation in (J2J) and (@J) is not important for us; it can be 
replaced by a more general one like 

p(a\knt-s\) 

with an integrable function p(r) decaying to zero as r — > oo. Since the exponential form seems to 
agree well with the Lagrangian measurements (see |37| for the Reynolds number around 100 and 
|41j for high Reynolds numbers) we will use it for the sake of simplicity. 
Set the rescaled velocity 

(6) U x (t,x) = X l - a U(X 2 ^t,Xx). 
Then U\(t, x) has the energy spectrum 

(nk) _( E (i-k®k\kr 2 )\k\i- 2 «, for \k\ e (to^e^x) 

However, we do not assume in this paper the full scale-invariance, namely, 
(8) U x (t, x) = U{t, x), for 4 = 0, £ = oo 

where = means the identity of the distributions. Instead, we assume the weaker assumption of the 
4th order scale invariance, i.e. that up to the 4th moments of the velocity field can be estimated 
in term of the energy spectrum as in the case of Gaussian fields. 

The viscous and integral scales £\,£o can be related to each other via the Reynolds number Re 

as 

£o _JL_ 
— ~ Re 4 -2<* 

£\ 

by using the positivity of kinetic energy dissipation of fluid in the limit Re — > oo. The correlation 
time a™ 1 ^! -2 ^ decreases as the wave number k increases. The spatial Hurst exponent of the velocity 
equals a—1 in the inertial range (£±, £$). It should be noted that because of the temporal stationarity 
of the Lagrangian field u(t, x + x^°'(t)) JS], [32] U(t,x) has the same one-time statistics as the 
Eulerian velocity u(t,x); in particular they share the same energy spectrum, but their multiple- 
time statistics are usually different. We could work with the modified von Karman spectrum but it 
is irrelevant for our purpose since we are concerned with transport in the inertial-convective range. 
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It is convenient to express the coefficients Eq, a in terms of Uq, the root mean square longitudinal 
velocity increment over the integral length £q, as 

(9) E Q « C a U^l~ 2a , a^CQlf- l U^ asV^^oo 
with dimensionless constants Co and 

_ (4^/ 2 2 2 »- 3 (2a-2)r(q + rf/2) 

(10) Ca ~ (d-l)r(2-a) 

where T(r) is the Gamma function. 

Assuming that the lifetime (i.e. correlation time r(k) = a -1 |/c| -2 ^) of eddy of size is same 
as its turnover time one gets the relation 

(11) a + 2/5 = 2. 

Assuming that the energy flux given by 8/^ fo)\k\/r(k) is constant across the scales in the inertial 
range one gets the relation 

(12) a -(3 = I. 

The values of parameters satisfying both eq. (|llj) and 1)12(1 correspond to the Kolmogorov spectrum 
with a = 4/3, f3 = 1/3. For the Kolmogorov spectrum, one has the expression, by estimating e by 

(13) E ^C a e 2/3 , a^c e 1/3 . 

Writing x(t) = — x^°\t) and adding the molecular diffusivity k we have the following Ito's 

stochastic equation for the pair separation x(t) 

dx(t) = [u(t,x(°\t)+x(t))-u(t,x(°\t))]dt + ^dw(t) 
= U(t, x(t))dt + ^dw(t) 

where w(t) is the standard Brownian motion in R rf . It is also useful to consider the associated 
backward stochastic flow which is the solution of the backward stochastic differential equation 

(14) d& a (x) = -U(s,& s (x))ds + y/Kdw(t), 0<s<t 

(15) **(*) = x. 

Denote by M the expectation with respect to the molecular diffusion and consider the scalar field 
T(t,x) 

(16) T(t,x) = M[T o ($o0*0)] 
which satisfies the advection-diffusion equation 

(17) ^^- = U(t,x)-VT(t,x) + ^AT(t,x), T(0,x)=T (x). 
We interpret eq. (|17j) in the weak sense 

(18) (T(t,-),9)-(T ,0) = ^ [ (T(s,-),A9)ds- f (T(s,-),V(s,-)-V9)ds 
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for any test function 9 £ C^°(M. d ), the space of smooth functions with compact supports. 

To study the small-scale behavior we introduce the following scaling limit. First we assume that 
the integral and viscous scales of the field U are £q = eL,l\ = e/K with L,K tending to oo in a 
way to be specified later. Then we re-scale the variables x — > ex, t — > e 2q t amounting to considering 
the re-scaled pair separation 

x *(t) = e^x^H). 
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The scaling parameter e will tend to zero, indicating that we are considering the emergent inertial 
range of scales l\ <C \x\ <C £q (since K,L — > oo) as a result of a large Reynolds number. We also 
set 

(19) k = e 2 ~ 2q K, with ii = «(e) 
After re-scaling, the advection-diffusion equation becomes 

(20) ^— = e^U^H, ex) ■ VT e + — AT 6 . 

\JTi Z 

We take the initial data T e (0,x) = T (x) G L°°(R d ) n L 2 (M d ). Let 

= e L - <x U{e lfi t,ex). 

As before (cf. ([7)l) the energy spectrum of the rescaled field V is given by 

Eb(I- A:® fcl/fcl" 2 )^! 1 " 20 , for \k\ G (L -1 , if) 



fc^(a,fc)- , 0> elge 
We rewrite eq. (|20[) in terms of V as 

(21) ^— =e 2q+a - 2 V(e 2{q - f3 h,x) ■ VT £ + ^AT £ , T e (0, x) = T (x). 

A simple, nontrivial scaling limit is the white-noise limit when 

(22) q < (3 
and 

(23) q = 2-a- P 

resulting from equating 2q + a — 2 and q — (3. Inequality (|2*2*|) and ((2*3*|) then gives the condition 

(24) a + 2(3 > 2. 
Note that for 

(25) a + (3 < 2 

and thus q > we have a short-time limit; otherwise, it is a long time (but small spatial scale) 
limit. 

The paper is organized as follows. In Section 2 we state the main results and discuss their 
implications. In Section 3 we discuss the meaning of solutions for the colored-noise and white-noise 
models and prove the uniqueness for the latter. In Section 4, we prove Theorem ^ we prove the 
tightness of the measures in Section 4.1 and, in Section 4.2, identify the limiting measure by the 
martingale formulation. In Section 5, we prove Theorem |2I The method of proof is the same as 
that in |14j (see also |S]). We refer the reader to [53] for the full exposition of the perturbed test 
function method used here. We note that the method of |23| requires sub-Gaussian behavior and 
spatial regularity of the velocity field and is not applicable here. 

2. Main Theorems and Interpretation 

Let us begin by briefly recalling the Kraichnan model. The model has a white-noise-in-time 
incompressible velocity field which can be described as the time derivative of a zero mean, isotropic 
Brownian vector field Bt with the two-time structure function 

E[B t (x) - B t (y)] [B s (x) - B s {y)} 

(26) = min (i,s) J 2[1 - cos (k ■ (x - y))]a _1 «S 1,(77+ 1, k^k^dk, 77 G (0, 1) 



with 



£ L {j] + 1, k) = lim £k,l{v + 1, k). 

if— >oo 



In this paper, we interpret the corresponding advection-diffusion equation for the Kraichnan model 
in the sense of Stratonovich's integral 



(27) dT t (x) = [VT t (x)Y o [dB t (x) - dB t (0)} + yAT t (i) dt, k > 0, T(0, x) = T (x) 



which can be rewritten as an Ito's SDE 

(28) dT t = f y A + ~sj T t dt + V2 a - 1/2 VT t ■ dW^ 
where W^(x) is the Brownian vector field with the spatial covariance 

(29) f( 1 \x,y)= J[exp(ik-x)-l][exp(-ik-y)-l]£ L (r ] + l,k)\k\ 1 ~ d dk, r) = a + /3-l 



and the operator B is given by 



(30) B4>{x) = £l1]W)|^g, <t>eCT<R*). 

hi 

We will discuss the meaning of solutions for the Kraichnan model and prove the uniqueness property 
in Section 3. The Kraichnan model for passive scalar has been widely studied to understand 
turbulent transport in the inertial range because of its tractability (see, e.g., [3B]> [TH] . [T2] . 
|30| . [27| . [2*U] . [TU] and the references therein). The tractability of this model lies in the Gaussian 
and white-noise nature of the velocity field. 

Theorem 1. Suppose a + 2/3 > 2. Let L < oo be fixed and let K = K(e) such that lim £ _»o K = oo. 
Let k = k(s) > such that lim £ ^o k = kq < oo. Let Tq G L°°(R d ) Pi L 2 (R d ). If, additionally, any 
one of the following conditions is satisfied: 

(i) a + 2/3 > 4; 

(ii) a + 2/3 = 4, lim £ ^ KeVlog-fT = 0; 

(iii) 3 < a + 2/3 < 4, lim £ ^ Ke 2 K 4 ~ a ~ 2/3 = 0; 

(iv) a + 2/3 = 3, lim e _>o ^e 2 K = lim £ ^o ^VlogK = 0; 

(v) 2 < a + 2/3 < 3, lim £ ^ K£ 2 if 4 ~ a ~ 2/3 = lim £ ^ eK 3 -°- 2 ^ = 

T/ien /or the exponent q given in \23\) the solution Tf of \21)) converges in distribution, as e — > ; 
in i/ie space D([0,oo);L™{R d )r\L 2 w (R d )) to the scalar field T t for pair dispersion in the Kraichnan 
model. The limiting Kraichnan model has the spatial covariance given by Here D([0, oo); L™* {R d )n\ 
L 2 u (R d )) is the space of L°°(R d ) n L 2 (R d ) -valued right continuous processes with left limits endowed 
with the Skorohod metric 3 and (R d ) (X^(lR d ) ) is the standard space L°°(R d ) (L 2 (R d ) ) endowed 
with the weak* (weak) topology. 

Remark 1. In addition to the assumptions stated in the Introduction and in the theorem, we use 
in the proof of Theorems Q] the assumption 

(31) sup/ \V t £ (x)\ dx = o ( -) , e^0, V0 < M < oo 

t<t J\x\<M \ £ J 

with a random constant possessing finite moments where 

1 s 
V t £ (x) = - E t£ - 2 V(-,x)ds. 
£ Jt £ 



For Gaussian velocity fields one has 



M a sup 

\x\<M 
t<t 



< CL a+2 ^ 2 log 



where the random constant C has a Gaussian-like tail by Chernoff 's bound. Condition h'Jl\) allows 
certain degree of intermittency in the velocity field. 

Note that, in Theorem^ when kq > and 2 < a + 2(3 < 3, lim £ _>o ke 2 K^~ a ~ 2 ^ = implies 
lim £ ^o sK 3 ~ a ~ 2 P = 0. Also, a + 2/3 < 3 contains the regime a + /3 < 2 in which the limiting 
Brownian velocity field is spatially Holder continuous and has a Hurst exponent rj = a + (3 — 1 G 
(1/2, 1), i.e. the limiting velocity field has a persistent spatial correlation. 

If we let L — > oo in the Kraichnan model, we see that it gives rise to a Brownian velocity field 
Bt with the structure function 

E[B t (x) - Bt(y)} [B s (x) - B s (y)} 

(32) = mm(t,s) J 2[1- cos (k-(x-y))]a- l B (a + f3,k)\k\ 1 - d dk 
where 

£{a + (3,k) = lim £ L (a + (3,k). 

L—>oo 

The spectral integral in ()32|) is convergent only for a + (3 < 2. The convergence of the integral in 
(|32l) means that the limiting Brownian velocity field Bt has spatially homogeneous increments. 

We can prove the convergence to the Kraichnan model with velocity field Bt in the simultaneous 
limit of e — > 0, K, L —* oo if additional conditions are satisfied: 

Theorem 2. Suppose a + (3 <2 and all the assumptions of Theorem^ (Thus, only regime (v) is 
relevant) except for the finiteness of L. Instead, let L = L(e) —* oo such that 

(33) limL 2 ( a+2 ^- 2 ) e = 0. 

Then the same convergence holds as in Theorem^ The limiting Brownian velocity field Bt has the 
structure function given by h'Jty) . 

Remark 2. In addition to the assumptions of Theorem^ (cf. Remark 1) we use in the proof of 
Theorems the assumption 

(34) sup/ \V t e (x)\ 2 dx < CL 2( - a+2l3 -V-, e -» 0, L -» oo, V0 < M < oo 
with a random constant C possessing finite moments. For Gaussian velocity fields one has 

sup / \V t £ {x)\ 2 dx < CL 2{ - a+2 P- 2 ^ flog , e -»■ 0, L -> oo, V0 < M < oo. 

t<t J\ x \<M V e / 

One sees t/wti Condition \°J4\) is in some sense more tolerant of intermittency than is. 

Due to the divergence-free property of the velocity field, the pre-limit scalar field satisfies the 
energy identity ([HU, Chapt. Ill, Theorem 7.2) 

(35) J |T/(x)| 2 dx + k J \VTf\ 2 (x)dxds = J \T (x)\ 2 dx 

provided that T G L 2 (R d ). From Qgft we have the estimates 

\\T?\\l < \\T \\l / \\T*\\ 2 Hl ds<(t + -)\\T \\l *>0 
jo K 



where || • is the norm of the standard Sobolev space of square-integrable functions 

with square-integrable 1st derivative. Thus the law of T £ is naturally supported by the space 
of continuous L 2 (IR rf )-valued processes which are also in L 2 OC ([0, oo); ff 1 (!R a! )). Following [3] we 
consider the space 

Q = £>([0, oo); L 2 w (R d ) n L%. (R d )) D L 2 w>loc ([0, oo); H 1 ^)) 

where the subscripts w and loc denote the weak and the local topologies, respectively. 

In the case of k > 0, kq > the above observation and the tightness argument for Theorem ^ 
and [21 then imply the tightness of Tf in the space f2. We have the following corollary. 

Corollary 1. // k > and T G L°°(IR d ) n L 2 (R d ) then the convergence holds in the space in 
the following regimes: 

Case 1: Let L < oo be fixed and K — > oo as e — > 0. 

(i) a + 2(3 > 4; 

(ii) a + 2(3 = 4, lim e _> e 2 ^/\ogK = 0; 

(iii) 2 < a + 2/3 < 4, lim e ^ e 2 K 4 ~ a - 2 ^ = 0. 

Case 2: Suppose a + (3<2<a + 2(3 and L, K — > oo as e — ► suc/j i/iaf 

lime 2 ^ 4 -"- 2 ' 3 = lim L 2 ^- 2 ^ = 0. 

In particular, 

(36) ||r || 2 -limsupE[||Tf|| 2 ] 



= liminf/c / E[\\VT^\\l]ds 
Jo 

> k / E[||VT S |||] ds > 0, f > 0, un/ess T s = 0, < s < t 

where zs i/ie solution of the corresponding Kraichnan model. 

In the case of k > 0, ko = and To G L 2 n the limiting Kraichnan model conserves the 
L 2 -norm of Tj. The energy identity (|35j) then implies 

l|7?lll< l|r ||| = ||r t || 2 , Ve >0, Vt>0 

which in turn implies lim e _^o II^TIb = 1 1 ^Z - * 1 1 2 - Hence the weak sense of convergence in Theorem ^ 
and [21 can be strengthened to the strong L 2 convergence. 

Corollary 2. If k = and T G L°°(R d ) n L 2 (M d ) i/ien toe convergence holds in the space 
D([0, oo);L 2 (R d ) nL~,(R d )) in toe respective regimes listed in Theorem^ and\^ In particular, 

||r ||| - lim K[||Tf |||] = 0, t>0 a.e. 

We see that in the context of Corollary 1 there is positive dissipation (|36f) while there is none in 
the context of Corollary 2. The conditions of the limit theorems set a constraint for the presence 
of positive dissipation: On the observation scale e, if the molecular diffusion k is of order e 2 ~ 2q , 
then there is always positive dissipation no matter how slow l\ vanishes. On the other hand if 
k <C e 2 ~ 2q (i.e. Ko = 0) and the dissipation is positive, then 

o n< *\ A- a -2(3 

l x = 0(e , v= — 

6 — a — 2(3 

with v G (2, oo) in the regime a + (3<2<a + 2(3 (cf. (|4()|)). An open question is whether there is 
a positive dissipation as e, k —* with l\ = at the outset. If there is, then the Kraichnan model 
(|2T|) is unlikely to be the governing equation of the scaling limit (if exists). 
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In the case of k = 0, a still stronger sense of convergence holds since now eq. (|21[) is of first order 
and any locally bounded measurable function <j)(T £ ) of the scalar field satisfies the same equation 
(|18|) with k = 0. The same argument for the proof of Theorem ^ and |2] will then yield the following 
result. 



Theorem 3. Assume the conditions stated in Remarks 1 and 2. Let k = 0, Tq,(()(Tq) 6 L°°(R ) n 
L 2 (R d ) where <j> is a locally bounded measurable function from R to R. Then T^,(j)(Tf) converge 
in the space D([0, oo); L^, (R d ) D L 2 (R d )) to the corresponding Kraichnan model in the following 
regimes. 

Case 1: Let L < oo be fixed and K —* oo as e — * 0. 

(i) q + 2/3 > 3; 

(ii) a + 2/3 = 3, lim £ ^ e v / log^ = 0; 

(iii) 2 < a + 2/3 < 3, lim £ ^ £-ftT 3 ~ Q ~ 2/3 = 0. 

Case 2: Suppose a + (3<2<a + 2(3 and L, K — > oo as e — > suc/i f/iai 

lim ^ 3 - Q ~ 2/3 = lim L 2 ^- 2 ) £ = 0. 

Remark 3. The assertions of Theorems^\^\^ and Corollaries^[\^hold true for random as well 
as deterministic initial data. 

When the parameters are in the regime a + /3<2<a + 2/3, by taking the expectation in the 
Ito's equation with the Brownian velocity field B t one sees readily that the longitudinal relative 
diffusion coefficient is given by 

(37) ^ + i^.f«(x,x).^ « -C-l p E \x\^ 
2 a \x\ \x\ a ™ 

for k « 1, ?7 = 1 - q = a + ft - 1 € (1/2, 1) 



with 



r (1) (x,x) = lim r (1) (x,x) = C~I p Eq\x 



L 



2*7 



^ 2(a + /?-l),_ 2(a + /3-l) 

(1 + - A — r— ; -) 1 - - A — r— ; -z (8) x x 

a — 1 a — 1 



where C Q+j g is defined as in (jlOj) . except with a replaced by a + (3. The exponent g is related to 
the exponent p in the expression for the mean square pair separation as follows: 

(38) E|x| 2 (t) ~ a - p E p t p , p= l/q = 1/(2- a-/?) 

up to a dimensionless constant depending only on q + /3. Expressions (|38|) and Q37|) can be viewed as 
the generalization of Richardson's i 3 and 4/3-laws, respectively. In general, p € (2, 00), indicating 
super-ballistic (i.e. accelerating) motion as a result of a scale-dependent relative diffusivity. 

We now remark on the range of scales for which Theorem [21 is proved and Richardson's laws can 
be reasonably interpreted. Let e be the scale of dispersion. Then the limit theorem holds in the 
range 



(39) e < min 



2(a+2/3-2) 



X \ 5-2c*-4,(3 

To, 



3=5=3^ > lf «o = 
4-a-2/3 - t ^ n 
if^o>0 



In the usual situation with £0 = 0(1) the range of scales covered by the limit theorem has an upper 
limit of 

(40) q, with 7g{[J/3^' 2); Sf«°>o for « + 2/3>2>a + /3, 

which is limited to the low end of the inertial range depending on a, /3, Kq. It is not clear whether 
this is physical or a technical matter. Qualitatively similar restriction of Richardson's laws in 
synthetic flows has been observed in numerical calculation (cf. [1], [15]). 



If we stretch the validity of 1)38(1 and (|37[) by taking the limit a — > 4/3, /3 — > 1/3 from within the 
valid regime, the resulting exponents are p = 3, 2r/ = 4/3 in accordance with Richardson's laws. 
On the boundary a + 2/3 = 2 the scaling exponent g should be given by 

(41) q = (3=\- a /2 

which also coincides with the limiting value of (|23|) . With (|41|) and K, L ^ oo, the solution of (|21jl 
converges to that of the advection-diffusion equation with the molecular diffusivity kq = lim e ^o K 
and the time-stationary, spatially Holder continuous velocity field V whose two-time correlation 
function is 

E[V(t,x)®V(s,y)] 

= / [exp (ik ■ x) - l][exp (-ik ■ y) - 1] exp (-a\k\ 2 ~ a \t - s\)£(a, fc)|fc| 1-<i dfc, a G (1, 2) 

jR d 

which has the self-similar structure 



E 



V(X 2p t, Ax) ® V (X 2f3 s, Xy)] = X 2a ~ 2 E [V(t, x) ® V(s, y)] . 



In view of the 4th order scale-invariance property it is reasonable to postulate the temporal self- 
similarity on the mean-square relative dispersion as kq — > 

E\x(t)\ 2 = f(E ,a)t 1 ^, 

which has the same exponent as the limiting case of (|38l) as a + 2(3 — > 2, where the unknown 
function / satisfies the relation 

f(E , Xa)X- l 'P = f(X- 2 E , a), VA > 0. 

Dimensional analysis with then leads to the relation 

E\x(t)\ 2 = c R c- 1 /^El /m t 1 /e, 

where Cr is the generalized Richardson constant. For (3 = 1/3 the exponent p is 1/3 as predicted 
by Richardson's t 3 -law. However, since the limiting velocity field is non-white-in-time, the notion 
of relative diffusivity is not strictly well-defined. Therefore the temporal memory persists on small 
or intermediate time scales and the notion of relative diffusivity does not describe accurately the 
process of relative dispersion on the boundary a + 2(3 = 2. (cf., e.g., [U, ^H] and jSOl)- 
Let us consider the regime a + 2(3 < 2. The correct scaling is to set 

(42) 2q + a-2 = or q = l- a/2. 

Then the exponent 2(q — (3) of the temporal scaling in (|21|) is positive due to a + 2(3 < 2, meaning 
the time variable is slowed down as e — > 0. It is easy to see by a regular perturbation argument that 
the solution Tf converges in the sense described in Theorem^ to the solution Tf of the following 
equation 

^ = V(0, x) ■ Vf t + ^Af t , f = T G L°°(R d ) 

if kq > 0. If, however, kq = 0, the above equation probably have multiple solutions for a given initial 
condition. The relation (|42|) is consistent with the numerical simulation using two-dimensional 
frozen velocity fields with Kolmogorov-type spectrum |llj . 

Unlike the previous regime, for either a + 2(3 = 2 or a + 2(3 < 2 there is no restriction on the 
vanishing rate of l\ . 
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3. Formulation 



From the general theory of parabolic partial differential equations ^7j> f° r an y fixed k > 0, e > 0, 
there is a unique C 2+,? -solution T £ (x), < Vr/ < a — 1. But the solutions T 4 e may lose all the 
regularity as k — ► 0, e — * 0. So we consider the weak formulation of the equation: 

(43) (Tf,6)-(T O ,0) = ^ j\ T !,M)ds-- £ f^(TlV{^-)-Ve)ds 

for any test function G C£°(R d ), the space of smooth functions with compact support. On the 
other hand the energy identity (j3SJ) implies Tf G L 2 ([0, to]; H 1 ^)) if T G L 2 (R d ). Hence for 
L 2 initial data the prelimit measure P e is supported in the space L 2 ([0, to]; F 1 (M d )) and, by the 
tightness result (Section 4.1), the limiting measure P is supported in L 2 ([0, to]; H^(M. d )). 
As in (|14|) and (|16[) the solutions T t £ can be represented as 



(44) T £ = M[T ($*' e (x))] 
where $s' e (x) is the unique stochastic flow satisfying 

(45) d& s < e (x) = -~V{-^,$l' £ (x))ds + Vkdw(t), 0<s<t 

(46) $t £ {x) = x. 

In the case of k = 0, $ ' e (x),Vt, is almost surely a diffeomorphism of R rf and T t e = To($g £ (x)). 



Moreover, for any locally bounded measurable function (p : R — > R, (j)(Tf(x)) = ((f) o To)(<l?g e (x)). 
In view of the averaging in the representation (|44|) we have 

Proposition 1. 

II^THoo — ll^olloo ft-S. 

Clearly, Proposition 1 holds for the case of k = as well. 

For tightness as well as identification of the limit, the following infinitesimal operator A £ will 
play an important role. Let Vf = V(t/e 2 ,-). Let Tf be the cr-algebras generated by {Vg, s < t} 
and Ef the corresponding conditional expectation w.r.t. T\. Let M e be the space of measurable 
function adapted to {.Ff , Vi} such that sup t<to E|/(t)| < oo. We say /(•) G V(A £ ), the domain of 
A £ , and „4 e / = gif f,geM £ and for / 5 (t) = $ _1 [E?/(< + 5) - f(t)] we have 

supE|/ 5 (t)| < oo 

t,6 

limE\f(t)-g(t)\ = 0, Vt. 

For f(t) = <p((T t £ ,e))J'(t) = <j)'((T t £ ,6)),V<f) G C°°(R) we have the following expression from @$ 
and the chain rule 

(47) A £ f(t) = ^f'(t)(T t £ ,Ae)-- £ f'(t)(T t £ ,V!(e)) 
where 

(48) Vf(9) = V t £ -V9. 
A main property of A £ is that 

(49) /(*)-/ -4 e /(s)ds is a Ff -martingale, V/eP(/). 

■/ o 

Also, 

(50) E^/(t) - /(a) = f\ £ s A £ f(r)dr Vs < f a.s. 

JO 
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(see El). 

Likewise we formulate the solutions for the Kraichnan model (|28|) as the solutions to the cor- 
responding martingale problem: Find a measure P (of Tf) on the space -D([0, oo); L^* (M. d )) such 
that 



(51) f((T t ,0))- j {f'((T s ,9)) 



^.(T S ,A9) + -(T S ,B*( 
2 a 



1 



+ ~f"((T s ,9))(6,iC%e) \ds 



is a martingale w.r.t. the filtration of a cylindrical Wiener process, for each f G C°°(]R) 
where B* is the adjoint of B and 

(52) (0,K®0) = J J T s {x)T s (y)V9(x)-T^(x,y)-V9(y)dy 
with T^\x,y) given, respectively, by (j^Hjl and 

(53) f (1) (x,y) = /[exp(iA:-x)-l][exp(-iA:-y)-l]f(?? + l,A:)|A:| 1 ^d/i:, ?? = a + /? - 1 



for L < oo and L = oo. To identify the limit for the proof of convergence one needs the uniqueness 
of solution to the martingale problem (|51jl which can be easily obtained as follows. 

Taking expectation of (|51jl with f(r) = r n ,n € N we get for the n— point correlation function 

F%(xi,x 2 ,x 3 , ...,x n ) = E To [T t (xi)Tt{x 2 ) ■ ■■T t (x n )] 

the equation 

{Fl® n 9) - {Fl® n 9) 



ds 



+ J2 - ( F n, f(1) ( x i, x j) ■ Ofa) ■ ■ ■ V0( Xi ) ■ ■ ■ V9( Xj ) ■ ■ ■ 9{x n ) 



t<3 



which induces a weakly continuous (hence strongly continuous) sub-Markovian semigroup on L p (lR' ld ),| 
\/p G (l,oo). The sub-Markovianity property is inherited from the pre-limit process Tf. The gen- 
erator of the semigroup is given formally as 

n i n 

(54) £ n $(xi,--- ,x n ) = ^Y, A ^ + -Y, T{1) ^ x ^ :V ^ V ^> k >0 



with the spatial covariance tensor (xi,Xj) given by (|29|) and (JEEJ), respectively, for L < oo and 
L = oo. Note that the symmetric operator C n (|54jl is an essentially self-adjoint positive operator 
on C^°(M> N ), N = nd which then induces a unique symmetric Markov semigroup of contractions 
on L 2 (R N ). The essential self-adjointness is due to the sub-Lipschitz growth of the square-root of 
(xi, x 2 ) at large \x\\, \x2\ (hence no escape to infinity) 
By Theorem 1.4.1 of 0| this semigroup induces a sub-Markovian Co-semigroup on L P (R N ), 
p G [1, oo). The uniqueness holds for these semigroups in their respective space as well but we will 
not pursue it here. 

li 



4. Proof of Theorem [T] 

4.1. Tightness. In the sequel we will adopt the following notation 

f(t) = f((T £ ,9)), f'(t)^f'((T £ ,9)), f(t) = f ({If, 9)), V/eC°°(R). 

Namely, the prime stands for the differentiation w.r.t. the original argument (not t) of /, /' etc. 

A family of processes {T £ ,0 < e < 1} C D([0, oo); (M d )) is tight if and only if the family 
of processes {(T £ ,d) ,0 < e < 1} C L>([0, oo); L^* [R d )) is tight for all 9 E C c °°(M d ). We use the 
tightness criterion of 25 (Chap. 3, Theorem 4), namely, we will prove: Firstly, 

(55) lim limsupP{sup | (T t £ , 9)\>N} = 0, Vi < oo. 

N->oo £ _>Q t <t 

Secondly, for each / £ C°°(M) there is a sequence f 6 (t) G £>(.A e ) such that for each to < oo 
{A e f e (t), < £ < 1, < t < to} is uniformly integrable and 

(56) limP{sup|/ e (t) - f((T £ ,9))\ > 5} = 0, V<5 > 0. 

e-+o t <t 

Then it follows that the laws of {(T £ , 9) , < e < 1} are tight in the space of D([0, oo); (R d )) 
Condition l|55j) is satisfied as a result of Proposition 1. Let 

/f(t) = i f°Ef/'(t)<Tf,V s e (0)) ds 
£ Jt 

be the 1-st perturbation of /(t). We obtain 



(57) fm = £ -f'(t)(T t £ ,V £ (9)) 
with 

(58) Vf{9) = Vf-V9 

(59) = v(l,^j=^ J°°W t V!ds 

where V has the power spectrum £k,l(o + 2/3, /c) by the spectral representation 

(60) E £ t V s £ = f[e ix - k -lje-^^-^Vfidk), Vs>t. 



Note that while Vf loses differentiability as K — > oo, V^ £ is almost surely a C '^-function in the 
limit with 

0<V??<a + 2^-2 
and has uniformly bounded local W^'^-norm, p > 1. 

Proposition 2. 

lim supE|/f (t)| = 0, lim sup |/i(i)| = in probability 

Proof. By Proposition ^ we have 

(61) n\fl(t)\] < £ ll/ , ||oo||To|| 00 ||^|| 00 / E\V t e \dx 

a J\x\<M 

and 

(62) sup|/f(t)| < -H/'lloollTollool^llooSup / \V{\dx. 

t<t a t<t J\x\<M 
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By the temporal stationarity of Vf we can replace the term E|V5- £ (x)| in (jtjlj) by E|V"(0, x)\. By 
assumption (cf. (|31|) . Remark 1), we have the desired estimate. Proposition |2] now follows from 
(jnij), dS2J) and (E2). □ 

Set / £ (i) = f(t) — fl(t). A straightforward calculation yields 

-47! = - ^/"(*) <7T, A #> {t!M{0)) + g/'(t) (r t £ , Av t £ (0) 
+ -/"(*) <T t £ , V £ (0)> (T{,V*(9)) - -f(t) (r £ ,v t £ (v t £ (0)) 
-~f'(t) <T £ ,V £ (0)> 



and, hence 
(63) 



<T £ , A(9) + -/'(*) (r £ ,Vf (V t £ (0))) - -/"(t) <T £ , V^)) <T t £ , V £ (&) 
2 a \ la 



+ 



2^ 



/"(t) <T £ , A0) <T t £ , V t £ (0)) - /'(t) (r £ , AVf (5) 



= ^(t)+Al(t) + A £ 3(t)+A|(t) 

where Af,(t) and A§(i) are the 0(1) statistical coupling terms. 

For the tightness criterion stated in the beginnings of the section, it remains to show 

Proposition 3. {„4 £ / £ } are uniformly integrable and 

limsupE|i4|(t)| = 

e^O t <t 



Proof. We show that {^4|}, i = 1, 2, 3, 4 are uniformly integrable. To see this, we have the following 
estimates. 

\Am\ = f !/'(*) (T t £ ,M)\ < fll/'lUllTolUIIAeil! 



Thus A\ is uniformly integrable since it is uniformly bounded. 

1 



\m)\ 



a 



f'(t)(T?,Vf(Vf(9)) 
T, 



< — ll/'lltx I 'u| X 

a 





1/2 




1/2 


[ \V t e \ 2 dx 




f \VV t £ \ 2 dx 




J\x\<M 




J\x\<M 





Similarly, 



\m)\ 



< 



a 
C 



nt) {if, vm)(T t £ ,vm 



f"|| \\T II 2 



\V/\ 2 dx + 



\x\<M 
13 



\x\<M 



\V t £ \ 2 dx 



Thus A| and A| are uniformly integrable in view of the uniform boundedness of the 4-th moment 
of Vf , Vf and Wf as L < oo is fixed and K — > oo (the 4th order scale- invariance) . 



I A £ \ 



KS 



< 



2a 
~2a 



\f"(t) (Tf,A9) (T{,V?(9)) - f(t) (T{,AVt(e) 

1/2 



|2 



(64) 



uriiooiiTo 

V t £ \ 2 dx + 



V t £ \ 2 dx 



x\<M 



+ H/'HooMoo X 



i|<Af 



ivv; £ | 2 dx + 



\x\<M 



x\<M 



\AV t £ \ 2 dx 



1/2 



The most severe term in the above argument as a result of K — > oo is 

K£ 

2a 

whose second moment can be bounded as 



f(t)(T t £ ,AV £ (9) 



f>(t)(T{,AVI(e) 



KB , 



(L E 


~\AVf\ 2 


dx 


\J\x\<M 







1/2 



(65) 



< Ci — H/'lloollTol 
2a 



K 3-a-2{3^ for q + 2/3 < 3 
< C 2 ke x < ^log K, for a + 2/3 = 3 
1, for a + 2/3 > 3 



and, thus, vanishes in the limit by the assumptions of the theorem. The 4-th moment behaves the 
same way by the 4th order scale-invariance. Hence A\ is uniformly integrable. Clearly 



limsupE|^!(t)| = 0. 

t<to 



□ 



4.2. Identification of the limit. Once the tightness is established we can use another result in 
.25] (Chapter 3, Theorem 2) to identify the limit. Let A be a diffusion or jump diffusion operator 
such that there is a unique solution oj t in the space D([0, oo); L^* (W 1 )) such that 

(66) f(u t ) - I Af(u s )ds 

Jo 

is a martingale. We shall show that for each / £ C°°(M) there exists f £ 6 T>{A £ ) such that 
(67) 

(68) 
(69) 

(70) 



sup E\f(t)-f((T t £ ,9))\ < oo 

t<to,£ 



UmE\f%t)-f((T{,e))\ 

e^O 



0, Vt < t 



sup E|^ e / £ (t)-^/((T t e ,0))| < oo 

t<to,£ 



limE\A £ f £ (t)-Af((T £ ,9))\ 



0, Vt < to- 



Then the aforementioned theorem implies that any tight processes {Tf,9) converge in law to the 
unique process generated by A. As before we adopt the notation f(t) = f((T£,6)). 
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For this purpose, we introduce the next perturbations /|,/|. Let 

(71) 4\cP) ee 

(72) 4\<p) = (4>,E[vt(ym) 

where the positive-definite operator K.} is defined as 



(73) 



= J 0(y)V</>(x)-TM(x,y)V<l)(y)dy 
(74) rW(x,y) = J[exp(ik-x) - 1] [exp (-tfc • y) - l]S KiL (a + p^^k^dk 

such that 



(75) 
(76) 

(cf. EH). 

It is easy to see that 
(77) 
(78) 



G 



(i) 

Ol,02 



4>(x)(f)(y)G^ e2 (x, y) dx dy 

d 2 



— ' dx l dyi 



{ {x)0 2 {y)Tfhx,y) 



4 1) (0)=E \(<f>, V e t {e)) (WW) 



4 1} (0 



where the operator B is given by 



d 2 4>(x) 
dx l dxi 



hi 



Define 



m) = -f(t) 

a 

/!(*) = \f{t) 



(:/v.v-(^)(:/; ; .v:(^}-A«(i?) 



ds 



Let 



G ( 2 ) fg ,|-Vr' 2) (x ) de 1 (x)de 2 ( y ) 
^oift^'V) - l^ L ij y x >y> dx i Q j 

9 u Kfe 2 )= 11 <f>(x)ct>(y)Gfl e2 (x,y)dxdy 



where the covariance function T^(x,y) = E Vf(x) (g) Vf(y) has the spectral density £k,l{ol + 
2/3, k). Let 

4 2) (0) = (0,Efe(V t ^))" 
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Noting that 

Ef[V s %x)®V £ (y)} 



(79) 



= J J [e ix - k - l}[e- iy - k ' - l] e ~ a \ k \ 2fi \ s -^ 2 e -a\k'\^\s-t\e 2 ve( dk ) g, Vf*(dk') 



+ y [e^ fe -l][ e -^' fe -l] [l 

we then have 
(80) 

and similarly 

(81) 

Proposition 4. 



-2a\k\ 2p \s-t\e- 



mo = ^rw 



£k,l(® + fc) cZA; 



Tf,Vf(0)) -4 2) (3f) 



/!(*) = l^/'W [(3?, vf (vf (*))) - 4 2) TO 



limsupE|/f(i)| = 0, limsupE|/f(t)| = 0. 



Proof. We have the bounds 

supE|/f(t)| < sup^HriUllTollLllV^IlL 



< C\e 2 



[ E\V t £ \ 2 (x)dx + [ \T^(x,x)\dx 

J\x\<M J\x\<M 



\x\<M 



supE|/f(t)| < sup-^H/'lloollTollc 
t<t t<t *a 



J\x\<M 



< C 2 e 2 K 2 ' a ' 2 ^; 
both of them tend to zero. 
We have 

1 



E\V t £ \ 2 (x) dx 



x\<M 



1/2 



E\VV t £ \ 2 (x)dx 



x\<M 



1/2 



□ 



A e fi{t) = -J\ t )[-{T £ ,Vtm{T £ ,Vl{0))+A^\T £ )\+Rl{t) 

A e m = -j\t)[-(rf,vt{vm))+4\Tt £ )]+Rm 



with 



m) 



fit) 



£ ^(T t £ ,A6)-^(T t £ ,V £ (9)) 



T!,Vt{6)) -4\T £ ) 



+nt)(T £ ,Vf(9) 



Re 2 



3?,Vf(VfW) 



(82) 



-f"(t) 



^ (Tf,AGfh t £ ) - -i (T £ ,Vt{Gf If) 



(2) 

where G e denotes the operator 



GfU = J GW{x,y)4>{y)dy, 
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and similarly 



he 2 



f"(t) 
+f(t) 
-fit) 



^<T/,A0)-^<T/,V^)> 



T t £ ,Vf(V!(e)))-A^(Tf) 



^,Avf(vr(^)))-^(/;-.V:(V / (V / (^)) 



Re 2 

{If, AE[Vf (V t e (0))] ) + ^ ( 7 )- . V/ ! E V; ( Vr (+))]) 



Proposition 5. 



limsupE|i2§(t)| =0, IimsupE|U|(t)| = 0. 



Proof. The argument is entirely analogous to that for Proposition 21 The most severe term without 
the prefactor k occurs in the expression for and can be bounded as 



< Cie||T ||oa ( I E\V t £ \ 2 dx 

\J\x\<M 



1/2 



(83) 



E 



\x\<M 



\vf\ 4 



E 



2t?£|4 



1 1/2 



dx + 



1/2 



E 



\x\<M 



ei4 



f/.X 



by assumption. The right side of the above tends to zero if either 

a + 2/3 > 3 

or 

(84) 
or 
(85) 



a + 2/3 = 3, limeydogif = 



a + 2/3 < 3, lim eif 3 -"- 2 ^ = 

e— >0 



is satisfied. The term involving e ( T t E , Vf (Cri T t e ) ) can be similarly estimated. 



The most severe term involving the prefactor k occurs in f?| and can be bounded as 



Re 2 E 



Tt^vmie)) 



< CRe 2 \\T< 



oo 



1/2 



E 



[x[<M 



3f/e|2 



(86) 



K£ 2 , 



for a + 2/3 > 4 
K£ 2 y/TogK, for a + 2/5 = 4 
- e 2^4- a -2/3 5 f or a + 2/5 < 4 



the right side of which tends to zero if either 

a + 2/5 > 4 

or 



or 
(87) 



a + 2/3 = 4, lim /ce 2 \/log if = 
3 < a + 2/3 < 4, lim Ke 2 if 4 ~ a ™ 2/3 = 
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or 

2 < a + 2/3 < 3, lim Ke 2 K 4 - a ~ 2(3 = lim eK 3 ' a ~ 2f3 = 0. 
Note that for a + 2(3 < 2 the condition (jHH) or (JSHJ) implies that 

hme 2 K i - a - 2 ' 3 = 0. 



Set 

W{t) = A%(t)-Rl{t)-Rm- 
It follows from Propositions 3 and 5 that 



□ 



limsupE|i? £ (t)| = 0. 



t<t 



Recall that 



M t %0) = f{t)- f A e f{s)ds 
Jo 



f(t)-fm-m)-m) 



'is 



f'(t) {Tf,A6) ds 



* 1 

o a L 



f"( S )4\T!) + f(s)A^(T!) 



ds 



R £ (s)ds 



is a martingale. Now that (|(i7|) - (|7()|) are satisfied we can identify the limiting martingale to be 



M t {9) = f(t) 



f'(s) 



Kq 

2 



(T S ,A9) + -A^(T S 



1 



where 



A? (</>)= lim 4 

K — >oo 



(1)/ 



a 

4 1} (0)= lim A 3 



K—>oo 



(cf. dj, O). 

Since (T t £ ,0) is uniformly bounded 

1(2?, 0>| < ||t o ||oo||0|Ii 

we have the convergence of the second moment 

limE{(Tf^) 2 }=E{(T t ^) 2 } 

Use f(r) =r and r 2 in (jKHjl 

M t (1) (0) = (T t ^)- 



— ■ Ci 



is a martingale with the quadratic variation 

2 



t a J 



A^(T s )ds = - 

a Jo 



MW(0),MW(0) 
where /CjP is a positive-definite operator given formally as 

(89) 4?0 = / %)vr t (x) -rw^.^vr^y)^ 

(cf. EH!). Therefore, 



(i) 



/c$!W s 
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where W s is a cylindrical Wiener process (i.e. dWt(x) is a space-time white noise field) and y KL Ta 
is the square-root of the positive-definite operator given in (|89jl. From (|7T|) and (|7H|) we see that 
the limiting process Tt is the distributional solution to the martingale problem (|51j) of the Ito's 
equation 

dT t = A + ^-Bj Ttdt + ^J 2^i£g dWt 

= (^A + ^T t dt + V2a- l / 2 VT t -dW t (1) 

where the operator B is given by (|3Uj) and W* is the Brownian vector field with the spatial 
covariance T^(x,y). 



5. Proof of Theorem [2] 

As we let L —* oo along with e — » the proof of the uniform integrability of A £ [f(t) — ff(t)] (the 
first part of Proposition 3) breaks down. In this case, we work with the perturbed test function 



nt) = f(t)-r l (t)+fi(t)+fi(t). 



Proposition 6. 



(90) lim supE|/f (t)| = 0, lim sup |/f (t)| =0 in probability, Vj = 1, 2, 3. 

Proof. The argument for the case of ff(t) is the same as Proposition 2. For /f (t) and /f(i) we 
have the bounds 



supE|/f(i)| < sup^liriloollTollLllVeil^ 



t<t 



t<to 

< Cte 2 L 2{ - a+2 ^- A 



\x\<M 



E\V t £ \ 2 (x)dx + 



\x\<M 



\t( 2 \x, x) \ dx 



supE|/f(t)| < Bup-jll/'HoollTol 
t<t 2a 



t<t 



|V0| 



\x\<M 



E\V t £ \ 2 (x)dx 



+ 



E\V t £ \ z {x)dx 



x\<M 



1/2 



E\VV t £ \ 2 (x)dx 



x\<M 



1/2 



< C 2 e 2 L 2 ^ a+2 ^- 4 

both of which vanish under the assumptions of the theorem. Here we have used the fact that 

/ E\V t £ \ 2 (x) dx = 0(L 2 ( Q+2 «- 4 ), L -» oo. 

J\x\<M 

As for estimating sup t<to \ fj(t)\,j = 2,3, we can use 



V t £ \ 2 (x) dx in place of 



E\V t £ \ 2 {x)dx 



x\<M 



'\x\<M 

in the above bounds and obtain by assumption (cf. (|34j) . Remark 2) the desired estimate which have 
a similar order of magnitude with an additional factor of 1/e and a random constant possessing 
finite moments. □ 
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We have 
(91) A £ f(t) 
with 
(92) 

and i?|(i), as before 
Proposition 7. 



(2?, A0) - -J"{t)4\T!) - \f{t)4\Tt) + R\(t) + Rlit) + R%(t) 



f'{t)(Tf,AVm 



limsupE|iS(t)| =0, j = 1,2,3. 

e ^°t<t 



Proof. The proof is similar to that of Proposition 5 with the additional consideration due to L 
These additional terms can all be estimated by 



oo. 



de 



E 



x\<M 



V t £ (x)V t £ (x] 



dx < CoeL 2 ^+ 2 ^ 



□ 



which tends to zero under the assumptions of the theorem. 

For the tightness it remains to show 
Proposition 8. {A £ f £ } are uniformly integrable. 

Proof. We shall prove that each term in the expression (|91j) is uniformly integrable. 

The first three terms are clearly bounded under the assumption of a + (3 < 2. The last three 
terms can be estimated as in Proposition by 

desup / V t E {x)V{{x) 

t<t J\x\<M 

whose second moment behaves like £ 2 _L 4 (° !+2 ^~ 2 ) ; by the 4th order scale-invariance property, and 
tends to zero. □ 

Now we have all the estimates needed to identify the limit as in the proof of Theorem ^ 
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